Inelastic quantum tunneling through disordered potential barrier 
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The effect of inelastic scattering on quantum tunneling through a rectangular potential barrier, 
of length L, containing randomly distributed impurities, is considered. It is shown that, despite 
the fact that the inelastic transition probability Tlnoiastic is exponentially small, it may be greater 
than that of the elastic in the out of resonance conditions for some energy interval AE (the energy 
of the incident particle is far from the spectrum of the random Hamiltonian) . The asymptotic of 
the tunnelling transition probability for various impurity configurations are found and the energy 
interval AE where T^iastic ~ exp(— L/L*) < Tlneiastic is determined. 
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I. INTRODUCTION 

Due to advances in manufacturing techniques of semi- 
conductor devices, that make it possible to create 
small tunnel junctions, and the need to understand 
fundamental conduction processes in order to improve 
these devices, quantum-mechanical tunneling of electrons 
through single, double or multi-barriers disordered reso- 
nant tunneling structure, is in the focus of physical inter- 
est Refii^ andi^^. The most important peculiarity of 
these systems is that some impurity configurations may 
increase the tunneling probability dramatically. At the 
same time, the change of any single localized state will 
change the magnitude of the resonant transmission co- 
efficient or even destroy the conditions needed for the 
resonant tunneling. There is the general belief that res- 
onant tunneling via a single localized state (or a single 
quantum well) can dominate the whole tunneling process. 
Assuming this, the resonant tunneling through the whole 
system can approximately be treated using the resonant 
tunneling via a single localized state. Then, the Breit- 
Wigner formula is postulated to express the transmission 
coefficient as 



T{E) 
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{E-Eor + {TL/2 + TR/2)^- 



(1) 



Here Eq is the energy of the localized state and F^, Tj^ 
are leak rates of an electron from the localized state to left 
and right leads respectively, taken as Tl oc exp(— 2x/L£) 
and Tji oc ex'p{—2(L—x)/Li), where L is barrier length, x 
is the position of impurity center and Li the localization 
length. The standard frameworks for the description of 
transport through a mesoscopic system are the Landauer- 
Buttiker formula, methods of random- matrix theory and 
Kubo's formula Refi ^i^i^ . Here we apply transfer matrix 
techniques, at first suggested in Refi^, with some gener- 
alization for inclusion of inelastic processes. 

When the energy of incident particle is far from the 
spectrum of the random Hamiltonian, which describes 
the barrier with impurities, than we are not in the res- 
onant conditions for tunneling. The Breit-Wigner for- 
mula is not applicable any more and the transition prob- 
ability is exponentially small T{E) oc exp(— L/L*) with 



L* = h/ ^/2m(U — E) where U barrier hight and L bar- 
rier length. The situation changes when inelastic pro- 
cesses are taken into account which is relevant for low 
temperatures when electron-phonon collisions become 
sufficiently inelastic. The incident particle, by emis- 
sion or absorbtion of a phonon, may be trapped into 
the resonant level with subsequent resonant tunneling. 
The probability of the inelastic scattering is exponen- 
tially small so that total tunneling probability is also ex- 
ponentially small but it turns out that in some energy 
interval, /\E = — L^ol, the inelastic tunneling proba- 
bility exponentially large compared to the elastic one in 
the regimes out of resonance. The last statement may 
be summarized in the following sequence of inequalities 

^olaltic(£^ - Eo) 1 » r(£;)i„elastic > exp(-L/L*). 

The aim of this work is to find out the energy interval 
AE when T(S)ineiastic ?J exp(— L/L*) and investigate 
the energy dependence of tunneling probability in that 
interval for various configurations of impurity centers. 



II. FORMALISM 

In this section we consider the tunneling of the quan- 
tum particle with the energy E through the rectangular 
barrier with ingrained in it randomly distributed impu- 
rities. In the case of short-range point-like impurities, 
randomly arranged in the points Zj (j — 1,2, ... , N), the 
impurity potential has the form 

N 

U{x)^U + Kj2^i^-^j)' 0<Zj<L (2) 

so that electrons are described by the single-particle 
Hamiltonian 



and interact with the phonon bath maintained in the 
thermal equilibrium through Frolich-like electron-phonon 
Hamiltonian. By introducing new variables 
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E)/h^, K = 2mKir? and fc^ = 2mE/h^ we arrive at the 
Schrodinger equation 



d^ 
da;^ 



N 



kH = lC^6{x- Zj) 



(4) 



with subsequent solutions of superposition the incident 
and reflected wave at a; < and transmitted wave at 
a; > i in the form 

\l/(a;) = exp(ifca;) + r exp(— iA;a;), a; < 0, (5a) 



and subsequent new boundary conditions 

V'+(0)(fc/«; - i) + ip-{0){k/K + i) = 2k/K, (13a) 



^+{L){k'/K + i) + i)-{L){k'/K - i) = 0. 



(13b) 



Starting from now we can construct the scattering matrix 
S that "carries" solution through the region < a; < i 
from the right to the loft boundary of the barrier: 



*(0) = S-^{L), 



(14) 



*(a;) = texp(iA;'(a; - i)), x>L. 



(5b) 



where r - reflection and t - transmission coeflicients. As- 
suming the continuity of wave function and it derivative 
at the boundaries of the barrier we have following bound- 
ary conditions 



*(o)+™=2, m=ik' 



ik 



*(L) 



(6) 



In the above defined variables the the barrier transmit- 
tancy T can be presented as 



(7) 



and it depends on the energy E of incident particle, num- 
ber of the impurity centers N and the configuration of 
impurities through the phase volume = {zj}. At 
this point it is convenient to introduce new notation for 
description the particle dynamics. Instead of the single 
wave function \E'(a;) we introduce the vector- wave func- 
tion 



*(x) 



V'-(a;) 



(8) 



such that 

^{x) = 'iIj_{x)+'iIj+{x), *'(x) = K(V'-(a;)-V'+(a;)). (9) 

Above definitions imply that the functions ip-{x) a^^d 
il>+{x) satisfy the following equation 



da; 



(10) 



everywhere, except the points zj so that 



d*(a;) 
da; 



+ i*(a;) = ^bJ2H^ - (11) 



where 



K 
-K 



B = 



1 1 
-1 -1 



(12) 



S = G{zi)n^2{z^\z+)G{Z2 - zi) . 



■ ■ Tn-\,n{z^\z'Ij)G{L — zn) 



(15) 



where G{x) is the free propagator that carries solution 
from right to left into the region free from the scattering 
centers 



QU) = ( exp(Ka;) 
^ I exp(— Ka;) 



(16) 



and T is the matrix that propagates solution through the 
scattering center 



Tj,j+i{zj\z-f) = 



1 



1-/C 



*(z, - 0) = f,.,+,{zj\z+)^{z, + 0). 
Also we introduce following two vectors 



(17a) 



(17b) 



that make it possible to present the transition probability 
in the very compact form 



TN{E,rN) = J 



\{f,S9)\^ 



(19) 



where function W{E, E') represents the probability of 
the electron transition from the quantum state with the 
energy E into quantum state with the energy E' as the re- 
sult of electron-phonon interaction, which may be found 
from linearized kinetic equation 



W{E, E') = W+{E, E') + W-{E, E') 



(20a) 



w- 



--{E,E')= J \M,f\- 



1 - npiE') 6{k{E) - k{E') t q) 



npiE) 



v{E') 



(20b) 
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x[{N„ + 1)6{E - E' - huj) + N„S{E - E' + huj)]p- 

where " ± " correspond to emission and absorption of the 
phonon, Mq - matrix element of the electron-phonon in- 
teraction, Nq and np - boson and fermion distribution 
functions and ^-functions presents conservation laws of 
the momentum and the energy. Representation Eq. l(T^ 
is very useful because all information about particle dy- 
namics in the barrier is included in S matrix, at the same 
time all information about boundary conditions lies in 
the vectors / and and inelastic processes included in 
the function W{E,E'). 

Caring out integration in Eq. (|20bl) . taking into ac- 
count only valuable exponential factors and comparing 
W to the elastic tunneling probability in the out of 
the resonance conditions oc exp(— L/L*) wc can find 
the energy interval Ai?, where inequality T(i?)ineiastic ^ 
exp(— L/L*) is satisfied, which equals to 



AE - Tin 



exp(-L/L*)-cxp(-|So|/r) 



(21) 



where T is temperature of the system and also assumed 
L>L*{\Eo\/T). 

Analyzing properties of S matrix it is possible to show 
that maximum value for the tunneling probability oc- 
curs for some particular impurity configurations (reso- 
nant configurations). The conditions needed for max{T} 



gives the probability of the resonant configuration real- 
ization Pj^{E) for the given energy of the incident particle 



"pres 

PNiE) = 

i N 



(25) 



Averaged over all possible configuration, with the fixed 
number of impurities, the resonant transmittancy has the 
form 



Tn{E)^-^ f TN{E,TN)dTN 



(26) 



and average value of transmittancy {T{E)) for arbitrary 
number of impurities can be found by averaging TJv {E) 
with Poisson distribution function (under assumption of 
absence of correlation between impurities) 



{T{E)) = PN^NiE), pjv = {cLf exp(-cL)/7V!, 

(27) 



N=l 



where c is linear impurity concentration. 



XiE,TN) = Sn-S22 - 1, YiE,TN)^Si2 + S2i ^ 1 

(22) 

define this resonant configuration with the statistical 
weight 



-pros 
N 



X(i3,zi...zjv)|<l 

r(i3,zi...i„)|<i 



dzidz2 ■ ■ -dz 



N- 



(23) 



It ratio to the full statistical weight of all possible con- 
figurations 



r 



N 



dzidz2 ■ ■ ■ dzN 



0<2l<...<ZjV<i 



-N 



Nl 



(24) 



III. TUNNEL TRANSPARENCY T{E) OF THE 
DISORDERED BARRIER 



A) One impurity scattering center A^ 1: In 

what follows we have introduced new parameters e = 
1 -I- /C/2k, C = Lk and t = ^2mL'^/{U ~ E). The f 
and 5* matrices in the main approximation equals to 



^^^1 1^1 ^ ^ I l-wr/2 2 + UT/2 



(28) 



eexp(£) exp(a;T(£ — Zi)) 



exp(— (£ — 2Zi)) exp{—uJT{C — Zi)) 



exp(£ — 2Zi) exp(wT(£ — Zi)) 2 exp(— £) exp(— a;T(£ — Zi)) 
The transmission probability has the form 



Ti{E, Tn=i) 
and subsequently 

{T^{E)) ^^j\^{E,Tj,=^)dZ^ -Z^^"" 



e^/T + 1 



1 + shih\C - 2Zi) + ujt{C - Zi)] + \e'^ exp(2/:) exp(2tjr(/: - Zi)) e^l'^ + e'^"/^ 

e^/^ + 1 



1 -I- exp 4/: 



r 



. l+LOT 
'2 + UJT 



exp -AC- 



.l+UJT 



2 + UJT J e^/T + ef'^/'^ ■ 



(29) 



(30) 



(31) 



B) Two impurity scattering centers = 2: con- Z2 > Z\ it is convenient to introduce following quan- 
sidering the case of two centers located in points Zi, 
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Coordinate z/L 



FIG. 1: In the graph is shown schematically the shape of the 
transparency profile as the function of single impurity position 
in the barrier for resonant elastic (solid line) and inelastic 
(dashed lines) scattering, with differen values of energy losses 



titles r] = Z2 — 2i - relative distance between centers 
and C. — {C — Z2) — Zi - asymmetry in the scattering 
centers configuration. Calculating elements of the scat- 
tering matrix S and then the X and Y functions form 
Eq. 1221) we have 

X{E, Tn=2) = [e^ - cxp(-2r7(l + ujt))] exp{C - Ziujt)- 

(32a) 

-[4 - exp(277(l + u;t))] cxp(-(L)) exp{{^C - Zi)ujt), 



Y{E,Tn=2) = 2sinh[C+(/:- 



Z2)wt] [eexp(77(l +ujt))+, 
(32b) 



+2exp(-77(l +tJT))] 

where rAr=2 = {??, C} < 77 < £ and \(^\ < C - rj. 
Analyzing the energy and configuration dependance of 
X{E,Tn=2) and Y{E,Tm=2) it is possible to determine 
several constraints, that define conditions for the max- 
imum tunnel probability and we list them below: a) 
e ~ exp(— 77(1 -I- ujt)) and esinh[C + [C — Z2)ujt] 1, 
b) for the energy interval exp(— £) < e < cxp(— £/2) 
we have 277(1 -I- ujt) ~ C and sinh[C + (£ — Z2)ujt] ^ 
exp(-/:/2), c) 2t]{1+ujt) - C and C+{C-Z2)ujt - 1 
for the energy e cxp(— £/2) and d) e ^ cxp(— £/2) in- 
dependently on configuration. By knowing 77 and C we 
can determine the probability for two impurity resonant 
configuration formation 



L'' J 0<ri<£. 



(33) 



and subsequently derive the resonant barrier transmit- 
tancy 



W) 



Ine 



'2+UJT 



2r{E) 



£[l + (2-e2/2)oxp(-£)] (^) 

(b) 



ln[eexp£/2] 2J^{E) 

2+ujT C(l+oxp(-£/2)) 



1-2UJT 



HE) 



l+2ujT J l+2cxp(-£) 

'■^^ exD (-C^^ 



(c) 



(d) 



(34) 

for (a) — (d) cases described above subsequently, with 
T{E) = (e^/^ + l)/(e^/^ + e^^*"/^). Also it is possible 
to estimate asymptotic behavior of the averaged barrier 
transparency for all possible impurity configurations 



(r2(i?))~— exp 



c n+ujT 

2 [l~UJT 



,E/T ^ 1 



(35) 

C) General case of N impurity centers arranged 




FIG. 2; One example of the shape for the transparency 
T2{E,Tm=2) profile as the function of impurity positions Zi 
and Z2- The distinctive maximum correspond to resonant 
configuration of type (c), by classification scheme for two- 
impurity tunneling, when the energy of the incident particle 
e ~ exp(— £/2). In the x and y axes is shown the ratio Z/C 
in pixel such that 10 pixels correspond to ratio Z/£ = 1 and 
in the z axis is shown absolute value of the tunnel probability. 

as < -El < . . . < Zjv < C: in this case detailed anal- 
ysis on the basis of above suggested formalism is very 
difficult, but rough estimation for the tunneling prob- 
ability {Tn{E)) may be found. It turns out that the 
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picture of tunneling under conditions of large A'^ ^ 1 
breaks into effectively two-centers problem. This is true 
if among all centers there are two, which separated as 
bZ = maxj^j+i — Zj} | In |e||. Such constraint comes 
from the fact that the 6*11 — ee^^ matrix element of 
the scattering matrix is exponentially dominant and if in 
some energy interval we can put ~ 1 than tunneling 
occur with maximum probability, otherwise it exponen- 
tially small. Subsequently, tunnelling via these centers 
dominate the whole tunneling process. The main contri- 
bution to the tunneling probability {Tjq{E)) in summa- 
tion formula Eq. (|27|l comes form optimal A'' such that 
N = A^opt ^ \J lnc| » 1. The resonant impurity con- 
figuration, with mentioned constraint for 8Z^ is in the 
order of - 2~^{N8Z - C)^'^/{N - 1)! and asymp- 



totic for the tunneling probability 

{Tn{E)) -exp(-2V/:|lnc| -£c- (ftw/T)V>C/|lnc|) . 

(36) 

To summarize, we have considered the tunneling of 
the quantum particle through disordered barrier in the 
regime of inelastic scattering when the tunneling proba- 
bility is greater, compared to that of elastic scattering, 
in the out of resonance conditions. The energy interval 
AE where such situation may be realized is found and 
set of the tunneling probability functions T{E) for 
various impurity configurations presented. 
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